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1 Introduction 

Two-dimensional Supersymmetric Quantum Mechanics [1] with supercharges of second 
order in derivatives was shown to be instrumental in building new integrable two-dimensional 
quantum and classical systems which do not allow the standard separation of variables [2] - 
[6]. All such system have, by construction, symmetry operators of fourth order in momenta 
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which, as a rule, are not reducible to lower order. Therefore, these models do not allow the 
conventional separation of variables (for definitions see [7]). 

The possible intertwining operator - supercharges - can be easily classified into four 
inequivalcnt classes depending on their metrics Qik [4]. The intertwining relations of second 
order can be represented in the form of a system of six nonlinear partial differential equations 
involving two partner potentials plus six functions defining the supercharge. It does not 
seem to be possible to solve this system in a general form. Till now, only some particular 
cases of constant metric were investigated [2] - [6]: elliptic (Laplacian) gn^ — Sik, Lorentz 
Qik — diag{l, —1), and degenerate — diag{l,0). The elliptic case leads to models with 
separation of variables and second order symmetry operators, but the other gave rise to a 
list of completely integrable models (with fourth order symmetry operators) which are not 
amenable to separation of variables. Some of these models were found previously by other 
methods [8], but some of them are new. We notice that although, in general, there is no 
direct relation between the integrability (with symmetry operators of order higher than two 
in momenta) and the solvability of the system, it was shown that some models with this kind 
of integrability allow partial [3] or even complete [6] solvability by means of SUSY methods. 

Only one class of intertwining supercharges with constant metric in the supercharges was 
not yet investigated: the case of a diagonal matrix gij, — diag{l, —a'^), a ^ 0, ±1. Just 
the study of this kind of intertwining relation is the objective of the present paper, the 
corresponding metric being naturally called a 'deformed hyperbolic' one. 

The paper is organized as follows. In Section 2 the solution of the intertwining relations 
with deformed hyperbolic metrics is performed, in such a way that the unknown functions 
satisfy unique nonlinear functional-differential equation. Since the general solution of this 
equation is very difficult to find, some suitable ansatze are investigated in Section 3, and 
particular solutions of the intertwining relations are obtained explicitly. In Section 4 the 
special case when variables are separated only in one of the partner Hamiltonians is con- 
sidered. A new SUSY-algorithm of separation of variables for the second partner potential 
is proposed for this particular class of models. In Section 5 we consider the solution of the 
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intertwining relations with deformed hyperbolic metric when one of the partner is chosen 
to be the isotropic harmonic oscillator. It is shown that its superpartner necessarily allows 
separation of variables too. Some final conclusions put an end to this paper. 



2 SUSY intertwining relations for supercharges with 
deformed hyperbohc metric 

Let us consider the SUSY intertwining relations 

HiQ+ = Q-Hi = H2Q- (1) 

between two two-dimensional partner Hamiltonians of Schrodinger type 

i/l,2 = -A + l/i,2(^) A = dj + di di = d/dXi f=(Xi,X2) (2) 

with second order supercharges of the form 

Q+ = Qikimdk + Ci{x)di + B{x) Q- = {Q+y. (3) 

These intertwining relations reaUze the isospectraUty (up to zero modes of Q^) of the su- 
perpartners Hi, H2 and the connection between their wave functions with the same values 
of energy: 

¥,'\x) = Q^¥^\x) ¥:\x) = Q-¥^\x) n = 0,1,2,... (4) 

Equations (1) are equivalent [2] to the following system of six nonlinear partial differential 
equations 

diCk + dkCi-iVi-V2)gik^0 (5) 
AC, + 2d,B + 2g,AV2 - - V2)Q = (6) 
AB + gikdidkV2 + CAV2 - {Vi - V2)B = (7) 
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where, for simplicity, the exphcit dependence of the functions on the variables has been 
eliminated. 

For the particular case gik — diag{l, —a^), a ^ 0, ±1, considered in this paper, the three 
equations in (5) take the form: 

29iCi = V 292C2 = -a'^v (8) 
diC2 + d2Ci = (9) 

where the notation v{x) = Vi{x) — V2{x) has been introduced. Prom Eq.(8) one can easily 
express the three functions Ci,2 and v in terms of a unique arbitrary function C{x) : 

Ci{x) = -\d2C{x) C2{x)^d,C{x) v{x)^-^did2C{x). (10) 

Then, due to (9), the arbitrary function C{x) must satisfy the second order wave equation: 

{df-^d^y{x)^o. (11) 

The new variables x± = xi± ax2 are obviously suitable to write the general solution of (11) 
in terms of two arbitrary functions C±{x±) as follows 

C(f) = Jc+{x+)dx+ + j C-{xJ)dx-. (12) 

Therefore, from (10) we get 

Ci(f) = --(C+(x+)-C_(a;_)) C2{x) = C+(x+)+C_(a;_) v{x) = --{C'^{x+)-C'_{x.)) 
a a 

(13) 

where the prime means derivative of the corresponding function with respect to its argument 
(we will also use below the notation d± = d/dx±). 

The two equations in (6) can be rewritten now in terms of the functions C+, C_ and their 
derivatives as follows : 

-^-^{Cl-C'^) + 2{d+ + d.){B + V2)-^{C'^-C'_){C+-C.) = (14) 

l±^(C^ + C'^)+2(9+-9_)(5-aV2) + ^(C;-C:)(C+ + C_) = 0. (15) 



Simple linear combinations of (14) and (15) lead to the system: 



ad+(2B + (1 - a'')V2) = -f (C; - Ci)C_ - (1 + a^)C'l - a(l + a^)d_V2 
ad42B + (1 - a'')V2) = |(C; - C'__)C+ + (1 + a^)C'; - a(l + a'')d+V2. 

The consistency condition for these two equations is: 



(16) 
(17) 



a^l + a^)V2 - {Cl + Cl) - a{l + a^){C' - C'_) 



0. 



(18) 



Therefore, the partner Hamiltonians can be expressed in terms of the four arbitrary functions 
C+(x+),C_(x_),Fi(xi) andF2(x2): 



Hi ^2 



-A + Vi^x) A = dl 
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+Fi(xi)+F2(x2). 



(1 + a^){dl + 9!) + 2(1 - a^)d+d. 



(19) 



(20) 



Prom (16), (17) and (20), one obtains the expression for B : 

B{x) = -lc+(x+)C_(x_)-i^(C;(x+)-Ci(x_))- 
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iClix+) + C!(x_)) - (Fi(xi) - a'F2{x2)). 



(21) 



2a2(l + a2) 

Thus, all these algebraic manipulations transformed the initial problem of solving the 
system of differential equations (5) - (7) to expressions (13), (20) and (21) in terms of arbi- 
trary functions C+(a;+), C_(a;_), Fi(a;i), ^2(2:2), which are restricted by the only remaining 
equation (7): 



d+{bC';{x+) + 2C+{x+)[Clix+) + Cl{x_) + 2F{x)]} = 
= d-{bC'l{x-) + 2C4x-)[Cl{x+) + C^_{x-) + 2F{x)]} 



(22) 



where function F and constant b were defined as 



F{x) = ^!i^t|!l(Fi(xi) - a'F2ix2)) b = -a'{l + a') 
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(23) 



The connection with previous works on the Lorentz metric case [2] - [6] can be estabhshed 
very easily. Indeed, Eq. (22) should be multiplied by (1 — a^), and the limit — > 1 should 
be taken. Then, the result given in Eq. (13) of [3] is straightforwardly obtained. 

Since Eq.(22) is the functional-differential equations for the functions C±^Fi^2-, which 
depend on their own arguments, there are no chances to solve this equation in the most 
general form. Nevertheless, in next section some particular solutions of (22) will be found 
by suitably chosen ansatze. In the following we are mainly interested in the models where 
at least one of the partner Hamiltonians is not amenable to separation of variables. 

3 Particular solutions of intertwining relations 

In order to obtain particular solutions of the functional-differential equation (22), some 
extra hypothesis must be done. We will consider the following three cases. 

3.1 Case in which C-{x-) = 

A great simplification is obtained if we choose one of the functions C± to be zero. Without 
loss of generality, let us consider C_(a;_) = 0. Then, one has to solve the second order 
equation for the functions C+{x+), F{x) : 

bC';{x+) + 2Cl{x+) + AC+{x+)F{x) = U.{x.) (24) 

where is an arbitrary function. Different possibilities appear here: 

1. U_{x_) = Canst 

From (24), this choice gives immediately that F[x) actually depends only on 
Therefore, from (23), Fi,2 (2^1,2) are at most linear functions of the corresponding ar- 
guments, and F{x) — c+a:+ -|- cq. This form of F{x) allows separation of variables 
for both Hamiltonians (19): to separate variables one has to rewrite (19) in terms of 
Zj^ — Xj^ — Xi -\- ax2 and — axi — X2- In the variables z± the kinetic part of the 
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Hamiltonian, in contrast to variables x± (see (19)), does not contain mixed terms, and 
the separation of variables for a linear function F becomes evident. 

2. C/_(x_) Const 

In this case, from Eq. (24), one can express F(x) as 

4F(f) - U+{x+)U_{x_) - L+(x+) (25) 

where the new functions are defined as: 



The special form (23) of the function F[x) means that 



(26) 



(dl - dV)F{x) = C/^(x+)C/_(x_) - U+{x+)U'L{x-) - L'[{x+) = (27) 
U'lix^) _ U'!!{x_) _ , 
C/+(x+) c/i(x_)-^ 

where 77 is an arbitrary real or purely imaginary constant. 

For ?7 = 0, the functions f/+(a;+), f/_(a;_) arc polynomials of first and second order, 
respectively. After their substitution into (27), and comparing with (26), one can check 
that both partner Hamiltonians allow separation, either in the variables x± or in the 
variables z±. 

For T] 0, the generic form of the functions U± is: 

U+{x+) = cr+ exp (770;+) + (5+ exp(— 770;+) (28) 
C/_(a;_) = (7_ exp(?7a;_) + 5_ exp(— ?7a;_) + 5 (29) 

where the constant coefficients a±, S± must assure the real chracter of U±{x±). Again, 
comparing (27) and (26), after a simple calculation one can obtain that L_|_ = br]^ 6 = 
0; ha+5+rf = 1/4, and 

Fi(xi) = 4q2^(x + 0?) (^+^"" 6xp(+^^^i) + ^+^- exp(-2?7Xi)^ + ki (30) 

-^2(3^2) = ~ AaHl + a^) (^+'^~ exp(+2a?7X2) + (5+(7- exp( -2077x2)^ + k2 (31) 



7 



where the constants ki are such that ki — a'^k2 = (1 — a^)rf/A:. The expressions (20) for 
the potentials Vi^2 include also the function C+(x+), which has to be found from (26), 
i.e. 

hC'[{x^) + 2Cl{x^) - WC+{,x+) = 0. (32) 

This equation can be integrated once: 

h{C'^{x+)f = - hrfl2f + C C = const. (33) 

For the arbitrary value of real constant C the function C+{x+) can be expressed in 
terms of elliptic functions (see [9]), and the corresponding potentials in (20) will be 
written in terms of such C+ and F from (30)- (31). Both partner Hamiltonians do not 
allow separation of variables. 

For the specific value C = 0, Eq. (33) becomes a couple of Ricatti equations [10] (let 
us remind that by definition h = — a^(l + a^)^ < 0): 

v^C;(x+) = ±{Cl{x^) - 67772). (34) 

The solutions are well known: 

C+{x+) = ^ ^^^^^^^P(^^^+/^) - ^exp(-i?7x+/y2) ^^^^ 
y/2 i/exp{i'r]x+/y/2) + iy exp{—ir]x+/y/2) 

where the constants u, i> have to keep C+(a;+) real. 

Choosing the lower sign in (35) one obtains from (20), (30) and (31) the following 
expressions for the partner potentials (up to a common constant term): 

V2{x) = Fi(a;i) + F2(a;2)-(l + aV/2 (36) 
= ]~ " a'((7+(7_e'''^i + 5+(5_e-'''^i) - (cr+^-e^"''^^ + 5+(7_e-'"''^^) 

V,{x) = V2{x)--C'^{x+) = V2{x)-r Cl{x+)-ril\ (37) 

Choosing the other sign in (35) is equivalent to the interchange Vi ^¥2- The constants 
in potential Vi can be tuned so that it grows at |a;| — > cxd, and the corresponding 
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Schrodinger equation will allow separation in terms of the variables xi,X2- The solutions 
of both one-dimensional equations can be represented in terms of Mathieu functions 
[9]. Therefore, the Schrodinger equation with potential Vi{x) in (37), which does 
not allow separation of variables, is isospectral (up to zero modes of supercharges) to 
the Schrodinger equation with potential V2 in (36). Since this last potential allows 
obviously separation of variables, one obtains a specific SUSY-separation of variables 
for potential (37). For further developments of this model, see Section 4 below. 

3.2 Two Ricatti equations 

The second ansatz which will give some solutions of the intertwining relations, i.e. solutions 
of (22), is: 

C±{x±) = cC|(x±) + d (38) 

with arbitrary non-zero real constants c, d. The general solutions of these Ricatti equations 
are: 

C± = d - — f± = cr± exp(7a;±) + S± exp(-7a;±) 7 = V-cd (39) 
J± 

where we can have cd > or cd < 0. 

Substituting relations (38) into (22), one can check that this last equation takes a simple 
form if ca(l -|- a^) = ±1, indeed : 

d+(^C+(x+)(^F(x) - d/c^^ = d_(^C_{x_)(^F(x) - d/c^ j (40) 

with F{x) defined in (23). This partial differential equation for F can be solved in general 
form (see [5]) by introducing the new variables t± = J dx±/C±{x±) : 

dt^C+CiF - d/c)) = dt_{C+C.{F - d/c)). (41) 

Therefore, its general solution is expressed in terms of an arbitrary function M of the com- 
bination (t_)_ + tJ). Hence, we have: 

F(f)-d/c = M(/^ + / ^)/C+(x+)C_(x_) (42) 
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where both functions C±{x±) are given by (39). Then, Eq. (42) reads: 

Fix) - d/c = lMQln(/+/_))/;/i//+/_ ^ [/(/+/_)/;/:. (43) 

Now one has to remind the special dependence of F{x) on x given in (23), i.e. did2F{x) = 0, 
which means that U"{f^f_) = 0. FinaUy, we obtain the following solutions for Fi 2 : 



^2(^2) = ^ (^a+5_e2«^^2 + (7_5+e-2"T^2) + ^ (^a^ 5^ e^"^^^ + a^^^e-^^^^^) + C2 



where ci — 0^02 — cd{l — a^). If we choose c = +l/a(l + a^), the partner potentials can be 
written in an explicit form 

V — Sda^S^ 
^ ~ a((7_eT^- - 5_e-'^=^-)2 ^ ^ 

' " ~ a{a+e^^+ - S+e-^^+f ^ ^^^^ 

up to a common additive constant. The alternative choice c = — l/a(l + a^) corresponds to 
the interchange Vi V2. 

We can observe that by choosing one of the following constants (T+, (T_} to be 
zero, we will obtain that only one of the partner Hamiltonians Hi, H2 admits separation 
of variables. For example, (7+ = allows to separate variables in H2 (but not in Hi)] 
the two one-dimensional potentials, that appear after separation of variables, are exactly 
solvable Morse potentials, and their eigenfunctions are known analytically. The separation 
of variables in H2 gives a chance to investigate a whole variety of bound states for the partner 
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Hamiltonian Hi too. In this sense one can speak of a specific kind of " SUSY-separation of 
variables" in Hi (for other types of SUSY-separation of variables see also [3, 6] and Section 
4 below). 

3.3 Case in which C'^{x±) = cC^{x±) 

This ansatz is a hmit case of the previous one, (38) with d — 0. The corresponding solutions 
are 

C±(a;±) = -l/cx±. (46) 
A calculation similar to the one carried out in the previous Subsection give us: 

F{x) = {x+x_)M{xl+x'i) M" = (47) 
-^1(2^1) = aixl + biXi F2{x2) — aixl + bio^xl (48) 
Vi,2{x) = ^^^±p- + ai{xi + xl) + bi{xi + xt). (49) 

The two isospectral Hamiltonians with potentials (49) do not admit separation of variables'^. 

4 SUSY-separation of variables 

As we noticed in the examples worked out in the previous Section, there are some models 
for which one of the partner Hamiltonians allows separation of variables, but the other does 
not. This means that the solution of two-dimensional Schrodinger equation for the separable 
Hamiltonian (for example, H2) is reduced to the solution of two separate one-dimensional 
Schrodinger problems. Its wave functions are expressed as a bilinear combination of these 
one-dimensional wave functions with arbitrary constant coefficients. We can use the SUSY 
intertwining relations (4) to obtain the wave functions of Hi, one can use ^'^^ from the 
now known wave functions \t'^^\ acting with the supercharge Q~^. But in contrast to the one- 
dimensional situation, some additional eigenstates of Hi may exist: if they are annihilated 
'^Due to the coefficients of tfie attractive singular terms in (49), tliese Hamiltonians are symmetric oper- 
ators, but they have no self-adjoint extensions (see [11]). 
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by Q , there are no partners in the spectrum of H2. Therefore, our task now is to find all 
the normalizable wave functions of Hi, which are simultaneously the zero modes of Q~. 

4.1 The algorithm 

To resolve this problem in the case of separation of variables in H2, we will use the following 
trick: though the variables are separated neither in Hi, nor in Q~ , we will consider a linear 
combination: 

Z = aHi + (3Q~ a, P = const. (50) 

Let us suppose now the existence of some constants a, (5, such that the variables in the 
operator Z are separated by some similarity transformation, and its normalizable eigenfunc- 
tions could be found. Then, one has to extract, among these eigenfunctions (by the direct 
action of Q~), those which are simultaneously the normalizable zero modes of Q~ . If this 
plan can be realized, we will reduce the spectral problem for the Hamiltonian Hi (which is 
not amenable to conventional separation of variables) to a couple of one- dimensional spectral 
problems. 

Keeping in mind [3] that sometimes the preliminary similarity transformations are helpful 
for separation of variables, we will transform the operator Z using a function exp{(/7(x)} that 
will be determined later on: 

Y = e-'^^^^Ze^^^^ ^{p- a)dl -{a + pa^)dl + 2{p - a){dnp)di - 
-2{a + (5a^){d2v)d2 - pC^dk + {P - a) ([dlcp) + (dicp)^'^ - 

-{a + Pa^) {{dl^) + {d2^f^ - pCkdkif + P{B - dkCk) + bVi. (51) 
To exclude the first order derivatives in Y, we have to impose two conditions: 

2{p - a){di^) = pCi -2{a + Pa^){d2^)=pC2 (52) 
which are satisfied if 

(a + Pa^)d2Ci = -(P - a)diC2. (53) 
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Comparing with (9), it is clear that the two constants introduced in (50) are related by 
2Q; = /3(l-a2). 

Substituting now diip,d2(p from (52) into (51), we get the expression: 
Taking into account that from (8) 9^(7^ = (1 — a^)(^i — V2)/2 and also (21), one obtains: 



I - a? I - a? 1 1 + 

B - d,C, + i-^K = 5 + '-^V, = -4C+C_ + ^(F, - F,] 



Now using (13), we can write (54) in the form: 

- ^(^■^ ^(-^ ^ ^ - ^) <-) 

Using the expressions (10) - (13) for Cj, the function (/^(x) can be written explicitly in 
terms of C±{x±) : 

'Pix) = ^(^"^2) (/ C+{x+)dx+ - jc.{x.)dx.y (56) 

Though for generic functions C±{x±) variables in (55) can not be separated, for some of 
ansatzc considered before this is just possible due to a specific choice of C±. In particular, 
the second ansatz of Subsection 3.1 with r/ 7^ 0, C_ = and C+ satisfying (34), eliminates all 
obstacles to separate variables for (55) in terms of Xi,X2- The same happens for the ansatz 
of the Subsection 3.2, if one of the constants cr+,5+ vanishes. In this case the functions 
C±{x±) satisfy also the Ricatti equation (38), and the potential V2 in (45) allows separation 
of variables. 

In both of these models, the similarity transformation above separates variables, and the 
eigenfunctions of the operator Y (and therefore, of Z) can be built from the eigenfunctions 
of the corresponding pair of one-dimensional problems. 

Now, we will briefly consider the particular case of the model with partner potentials 
(44) - (45) for (7+ = 0. Then, the physical system is described as 

H2 = h^{Xi) + h2{x2) hi = -dl + F^{xi) h2 = -dl + F2{X2) (57) 
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hiipi{xi) = eiV'i(xi) h2'ip2{x2) = e2'ip2{x2) (58) 

where the two one-dimensional Hamihonians hi and /i2 are, up to additive constants Ci,C2, 
Morse potentials with well known bound states, if (T_(5_ > 0. For simplicity, we will consider 
(7- — 5- = (7 and 0? < 0. Substituting (39) into (13) and (21) one obtains the expression for 
the supercharge : 



2u , , y-M^ + a')^ , (l-a2)ci 



Q+ = -hi + a'h2 + Cidi + ^ -C_ + 



a? a 



The action of this operator on the wave functions ^{x) — ipi{xi)ijj2{x2) of H2 will give the 
wave functions (if normalizable) of the partner Hi : 

/ o fl - a'^)d\ J-dail + a^) 

g+*(f) = (a^e2-ei+ ^ ^ )^(f) - ^ ^ '-{di - ad2)^{x) + 

\ a J a 

C ix ) f J-dail + a^)\ 
+ ':±i}^(d^^ad2-^ ^ -]^{x). (59) 



a 

The one-dimensional wave functions "01 (2^1) ^^'^ '^2(2^2) are known explicitly (see, for ex- 
ample [3]) in terms of hypergeometric functions. One can straightforwardly check that the 
singularity of C_(x_) at x_ = in the last term of (59) can not be compensated. Therefore, 
the bound states '^{x) of H2 have no normalizable partner states in the spectrum of Hi. 
According to the discussion carried out in the the first part of this Section, the additional 
bound states of Hi could be found among normalizable zero modes of Q~ . In the considered 
model, the operator (55) takes a simple form (up to a constant factor and a constant additive 
term), which allows the separation of variables: 

+ (60) 

Then, the eigenfunctions ^z(^) of the operator Z, defined in (50), are: 

= e'^(^)*y(f) = e'^(^) • Vi(a;i)^/^2(a;2) (61) 

where '0i(2^i)) '^2(2^2) are arbitrary eigenfunctions of hi,h2, respectively, and the similarity 
transformation is performed by 

sinh(7a;_) 
14 



The singularity of (62) at = is responsible of the non-normahzability of all eigenfunctions 
^z{x), including possible zero modes of Q~. Therefore, the Hamiltonian Hi has no bound 
states at all, despite the presence of two Morse potentials in (44). 



5 Two-dimensional harmonic oscillator 

In the previous successful attempts [2] - [6] to solve the two-dimensional intertwining 
relations of second order (1), the system of nonlinear differential equations (5) - (7) was 
solved step by step, and only at the very end some specific expressions for potentials were 
obtained. Up to now, it has been impossible to solve these equations straightforwardly, 
starting from the fixed form of one of the partner potentials, apart from some very simple 
examples with separation of variables in both Vi and V2. The soluition of the problem is not 
so evident in the case of supercharges with deformed hyperbolic metric. In this Section, we 
will take ^(x), as an isotropic harmonic oscillator: 



One can check that non-isotropic oscillator does not produce any solution of the intertwining 
relations. Substituting this potential in (18), we obtain 



where A is an arbitrary constant, and b was defined in (23). The variables are separated, 
and the functions C±{x±) satisfy the Ricatti equations: 



Then, equations (16) - (17) give the expression for the function B{x) in terms of C_|_, C_ : 




(63) 



dl [Cl + V^C\] = dl [C! - V^C'_] = 2A const 



(64) 



Cl{x+) + V^C'_^{x+) = A^x\ + A+X+ + a+ 
C!(x_) - ^^^C'_{xJ) = A^xi + + a_. 



(65) 
(66) 



/x(l + a^)+(x;(l-a2) 
2 



{x\ + xi) + 
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A-a2(/x(l-a2)+u;(l + a2)) 
+ x+X- + P+X+ + p-X- + const. 

After rather cumbersome but straightforward calculations, one can check that the last 
equation to be solved (7), is fulfilled only if, in addition to (65) - (66), one of the functions 
C+, C_ is linear, for example C_(a;_) ~ ax^. Then, we finally obtain that 

yj^) = yJ^) - -Ci(x+) + const. (67) 

a 

The Schrodinger equation with such potential Vi does not allow separation of variables in 
terms of x± due to the presence of a mixed term in the Laplacian, but it obviously allows 
separation in the variables z+ — x+ — Xi + ax2, — axi — X2, due to the absence of mixed 
term. Therefore, no non-trivial isospectral superpartners of the isotropic harmonic oscillator 
can be built in the framework of second order intertwining relations for superchrages with 
constant metric. 

6 Conclusions 

The investigation of intertwining relations (1) with deformed hyperbolic metric — 
diag{l, —a^), a ^ 0, ±1 in the supercharge operator (3), that has been carried out in this 
paper, completes the study of second order intertwining relations with constant matrix gik 
in two-dimensional Quantum Mechanics, which was started in [2] - [6]. In the present case 
also some particular classes of solutions for the partner potentials were found (see Section 3). 
Among them, there are pairs in which one of the partners allows separation of variables, but 
the second one does not. A specific procedure of SUSY-separation of variables is proposed for 
this case (see Section 4) . In a particular model, this new algorithm led to the conclusion that 
the system with the attractive potential (44) for (7+ = 0, (7_ = 5_ does not allow any bound 
states. Although the nonlinear equation (18) is not amenable, as a rule, to solution with 
given potential (e.g. for the Coulomb potential), sometimes this is possible. Thus in the 
case of one partner being harmonic oscillator, it is shown in Section 5 that the second one also 
allows separation of variables. All models considered in the paper are completely integrable. 
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i.e., nontrivial symmetry operators of fourth order in momenta {Ri — Q^Q~] R2 — Q~Q^) 
exist for the deformed hyperbohc metric also. 

A few additional ansatze could be considered in the same manner. In particular, one can 
check that the case C-{xJ) = Const 7^ (the natural generalization of the Subsection 3.1) 
leads to a particular sohition, obtained above within the ansatz in the Subsection 3.2. The 
partner of the constant potential V2 — Const also allows separation of variables, similar to 
the case considered in Section 5. In this paper we were only interested in real potentials. 
Some models with complex potentials (see [12], [6]) can be considered if for example we allow 
purely imaginary a 7^ ±i or arbitrary values of u, i> in (35). 
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